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ABSTRACT 



We present two-dimensional (2D) gas kinematics and excitation of the inner 300 pc 
of the Seyfert galaxy ES0428-G14 at a sampling of 14 pc^, from near-infrared spec- 
troscopic observations at Rw6000 obtained with the Integral Field Unit of the Gemini 
Near-Infrared Spectrograph. From measurements of fluxes and profiles of the emission 
lines [Feii]A1.257/im, Pa/3, Il2A2.121/im and Br7, we construct 2D maps of line intensi- 
ties and ratios, radial velocities and velocity dispersions. Emission lines "tomography" 
is provided by velocity slices obtained across the line profiles, a unique capability of 
IFUs, which allows the mapping of not only of peak velocities but including also the 
wings. We compare these maps with a previously published high spatial resolution 
radio map and find a tight relation between the radio structure and the emission-line 
flux distributions and kinematics, revealing that the radio-jet plays a fundamental role 
not only in shaping the NLR but also in the imprint of its kinematics. Blueshifts of up 
to 400kms~^ and velocity dispersions of up to 150kms~^ are observed in association 
with the radio jet at position angle PA ~ 129°, which is also the PA of the photometric 
major axis of the galaxy. We conclude that the radio jet is launched at a small angle 
relative to the galactic plane, with the NW side slightly oriented toward us. This angle 
is small enough for the radio jet to shock and compress the gas in the plane of the 
galaxy, and for the nuclear continuum to ionize and heat it. The distinct kinematics 
and flux distributions observed for the different emission lines suggest different origins 
for their emission. The [Fe ii] shows the largest blueshifts and velocity dispersions and 
its flux distribution is concentrated along the jet, while the H2 shows the lowest veloc- 
ity dispersions and has additional flux contribution from regions beyond the jet. Both 
X-rays emitted by the active galactic nucleus and shocks produced by the radio jet 
can excite the H2 and [Fe 11] emission lines. We use the 2D velocity dispersion maps to 
estimate upper limits to the contribution of the radio jet to the excitation of [Fe 11] and 
H2 which may reach 90% for [Fell] and 80% for H2 in the jet region. The [Fell] /Pa/3 
emission-line ratios and the association of the [Fe 11] flux distribution and kinematics 
with the radio structure supports a stronger contribution of the radio jet to the [Fell] 
excitation than to that of H2. In the regions beyond the jet the observations favor 
X-ray excitation. 
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1 INTRODUCTION 

The Narrow-Line Region (NLR) of Seyfert galaxies is 
one of the best probes of the mechanisms in opera- 
tion in the surrounding of accreting supermassive black 
holes in galaxies. The excitation and dynamics of the 
inner NLR gas can reveal how radiation and mass 
outflows from the nucleus interact with circumnuclear 
gas. Until recently, the measurement of these proper- 
ties had to rely on either n arrow-band images (e.g. 

[ Wilson et.aPllQQal: ICaoetti etll 1996,: Schmitt fc KinnevI 

Il996f) or long-slit spectroscopy. Long-slit spectroscopic 
studies provide information on the origin and excita- 
tion o f the emission lines (e.g. IVeil leux. G oodrich fc Hilll 
Il997l: iRodrfguez-Ardila. Riffel fc Pas toriza '2 00j) as well 
as on the gas kinematics (e. g. i,Winge_eL_al 



19991: ICrenshaw fc Kraemed I2OO0I: iHutchings et alj 



Ixaiser et alJl2000HNelson et al.ll20oi 
the locations covered by the long-slit 

, iMulchaev. Wilson fc Tsveta~ 

iFerruit. Wilson fc MulchaevI 12 



1997 



199S; 



but are restricted to 



and 

have shown that 
obscuration can affect the optical morphology of the emit- 
ting gas region, a problem that can be alleviated by using 
infrared lines to map the NLR. Relevant emission-lines in 
the near-IR include [Fe ll]A1.257/im and 1.644/im, HI lines 
such as Pa [3 and Br 7 and molecular hydrogen lines such as 
H2A1.957/jm and H2A2.121^m, w hich can be used to map 
the gas kinematics and excitation dStorchi-Bergmann et all 

Frngilwinge oHMmI. 

Studies based on the above lines iForbes fc Ward| 
I1993I : iBlietz et al.l Il994) have revealed a correlation be- 
tween the [Fen] and the radio emission, indicating that 
shock excitation by the radio jets is a likely mechanism 
for production of the [Fell] emission in these objects, 
although other works have favoured phot o ionization 
by the nu clear source ISimpson et alJ Il996l : iDas et all 
20051 . l2006h . H2 emission is al so strong in AGNs (e .) 



Veilleux. Goodrich fc HilJ Il997l: iR.odriguez-Ardila et all 



20041: iRodriguez- Ardila. Riff el fc Pastorizal |2005^ and can 

be used to probe the distribution of warm molecular gas. 
To further progress in this area and investigate the origin 
of the near-IR emission lines it is necessary to spatially 
resolve the kinematics and excitation of the emitting gas. 

The multiplexing capability necessary to map, simulta- 
neously, the gaseous distribution, excitation and kinemat- 
ics at such small scales with sufficient spectral resolution 
is now provided by the integral field units (IFU) in oper- 
ation in large telescopes. In this work we use the IFU of 
the Gemini Near-Infrared Spectrograph which provides two 
dimensional (hereafter 2D) mapping in the near-IR (thus 
minimizing the effect of reddening) to study the circumnu- 
clear emitting gas of the Seyfert galaxy ES0 428-G14. We 
present 2D flux, radial velocity and velocity dispersion maps 
in four strong emission lines: [Fe 11] A1.257/im, Pa/3A1.282/im, 
H2A2.121^im and Br7A2.166^im. The R?5i6000 spectral reso- 
lution has allowed us also to obtain velocity slices across the 
emission-line profiles in order to better probe the velocity 
fields. 

We present in addition broad-band J and K images and 
a narrow-band image in the line [Fe n]A1.257/.im obtained 
with the CTIO Blanco Telescope. We explore the relation 
between reconstructed images in the near-IR emission lines. 



their kinematic maps, the narrow-band [Fe 11] image and pre- 
viously published optical [O iii]A5007 narrow-band and radio 
images, in order to investigate the excitation mechanism of 
the near-IR lines, in particular [Fell] and H2. 

ESQ 428 -G14 is an SO galaxy which was classified as 
Seyfert 2 by iBergvall. Johansson fc OlofssonI lll98d) . High 
resolution VLA maps show predominantly a two-sided and 
asymmetric radio jet oriented approximately along the 
galaxy line of nodes at position angle PA« 129°, in which 
the NW side is less extended and terminates in a bright 
hot spot at 0'.'75 from the nucleus, while to the SE the 
radio emission is fainter, the je t bends to N and its ex- 
tent is at least three times longer ([^rosta^g Wilson Il989l : 
IFalcke et al. ll996l:lFaicke. Wilson fc Simpson 1 9981. [Oiii] 
and Ha-|-[Nll] images obtained with the Hubble Space Tele- 
scope show extended emission well aligned with the ra- 
dio jet with stronger emission to the NW as observed in 
the radio. The ratio [O ni]/(Ha-|-[Nll] ) shows a bipolar 
structure with larger values to t he SE jFalcke et al. Ill996l : 
IFalcke. Wilson fc Simpson ||l998l). We ad opt a distance for 
this galaxv of 19 Mpc fFalcke et anil996l) for which 1" cor- 
responds to 92 pc at the galaxy. 

This paper is organized as follows: in Section 2 we de- 
scribe the observations and data reduction. In Section 3 we 
present the near-IR spectra and images, a comparison with 
previous optical narrow-band and radio images and the kine- 
matics of the gas. In Section 4, we discuss the results and 
their implications for the origin of the [Fe 11] and H2 emission 
and in Section 5 we present the conclusions of this work. 



2 OBSERVATIONS AND REDUCTIONS 
2.1 Integral Field Spectroscopy 

The spectroscopic data were obtain ed with the Gemin i Near 
Infra- Red Spectrograph (GNIRS) iElias et al.l Il998l) Inte- 
gral Field Unit (IFU) in December 2004 under the in- 
strument science verification program GS-2004B-SV-26, and 
comprise two sets of observations centred at A=2.21/.im and 
1.27/im. We have used the 111 1/mm grating with the Short 
Blue Camera (0'.'l5/pixel) which gives a resolving power of 
R=5900. The GNIRS IFU has a rectangular field of view, 
of approximately 3'.'2 x 4'.'8, divided into 21 slices. At the 
detector, the slices are divided along their length into 0'.'15 
square IFU elements. For our observations, the major axis of 
the IFU was oriented along the PA « 129° , which is the ori- 
entation of the radio jet. The observing procedure followed 
the standard Object-Sky-Sky-Object dither sequence, with 
off-source sky positions since the target is extended, and 
individual exposure times of 600s. One of the IFU slicing 
mirrors was damaged during assembly and presents a lower 
(20% of the nominal value) transmission region, about l'.'8 
in length, offset 0'.'5 from the center, therefore small spatial 
offsets in both directions were added between exposures to 
"fill in" the signal in this region. Telluric standard stars were 
observed immediately after the target. The basic observing 
information is shown in Table Q Conditions during the ob- 
servations were good for the night of Dec. 27 (clear skies, 
image quality (IQ) ~0'.'60), and patchy clouds (but stable 
guide counts) with similar IQ for the Dec. 28 dataset. 

The data reduction was accomplished using tasks con- 
tained in GNIRS package which is part of the GEMINI IRAF 
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Table 1. Observations 



Instrument / 


Date 




Total Exp. 


Telescope 


(UT) 


Filter 


time (s) 


GNIRS/Gemini 


2004/12/27 


K_G0503 


2400 


GNIRS/Gemini 


2004/12/28 


J_G0505 


1200 


OSIRIS/Blanco 


2000/02/15 


K 


80 






J 


80 






1.06 


270 






[Fen] 1.257 


1800 



package as well as generic IRAF tasks. Through the reduc- 
tion tasks we have performed trimming, flat-fielding, sky 
subtraction, wavelength and s-distortion calibrations. We 
have also removed the telluric bands and flux calibrated 
the frames using the star HR4023 (A2V spectral type) as a 
"relative" flux standard. The cosmic ray cleaning was done 
before the sky subtraction using the algorithm described by 
Ivan DokkumI (i2001) . In Fig.Qwe present the nuclear spectra 
and the spectra at C/.'SNW from the nucleus, which is the 
position of the peak of the radio continuum emission within 
an aperture of 0'.'45 diameter. 



2.2 Imaging 

Broad-band J and K images and narrow band images cen- 
tered on the emission line [Fell]A1.257 and adjacent contin- 
uum were obtained in February 2000 at the Cerro Tololo In- 
teramerican Observatory (CTIO) using the 4-m Blanco tele- 
scope. The observations were performed with the infrared 
imager/spectrometer OSIRIS using a 1024 x 1024 HgCdTe 
array at an angular scale of 0'.'l6pix~^. Details of the ob- 
servations are presented in Table Q The angular resolution 
was 0'.'8-0'.'9 (FWHM of standard star images). 

The observing procedure consisted of taking a sequence 
of 4-9 dithered images on source with sky images taken be- 
fore or after the object sequence. Individual sky and object 
images were first divided by the fiat-field. We then median 
combined the sky images and subtracted them from the cor- 
responding object images. The individual images in the same 
filter were then aligned using as many common point sources 
(stars and the nucleus) as possible and finally averaged. We 
have used the task imcoadd from the GEMINI package to 
perform this alignment/averaging procedure. 

Photometric calibration was performed for both broad 
and narrow band filters using standard stars. Correction for 
extinction was performed u sing the mean extinction coeffi- 
cients for February llFrogellll99ah . 

The narrow band [Fe ll] image was calibrated by multi- 
plying each filter transmission curve by the black body spec- 
trum corresponding to the associated standard star spectral 
type and then using this integrated fiux as the fiux of the 
standard star in the filter. After calibration, the emission line 
image was created by subtracting the aligned continuum im- 
age (Filter 1.06/im) from the continuum plus emission line 
image. 
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Figure 1. Top: Nuclear spectra for ES0 428-G14 in the J-band 
(left) and K-hand (right). Bottom: Spectra at 0''8NW from the 
nucleus, which is the position of the radio continuum emission 
peak. 



3 RESULTS 

In this section we discuss the morphology of the narrow- 
band [Fell] image and the color map J~K constructed from 
the broad band images, comparing our data with narrow- 
band optical [O III] and radio images available from previous 
studies. We also present the 2D maps obtained from the IFU 
spectroscopy: emission-line distributions and ratios, velocity 
field and velocity dispersion maps. 



3.1 Images 

Narrow-band [Fell] image: In panel a of Fig.|21we present 
our [Fe ll] image which shows emission extended by ~ 3" to 
the SE and 1'.'5 to NW along PA« 129°. For comparison, 
we show also in Fig.|21 an HST narrow band [O ill] image 
(panel b) as well as a radio 2-cm image (panel c), both from 
iFalcke et all lll99efl . The [Fell] image shows a double "V"- 
shaped morphology, possibly indicative of a bi-cone. There 
is also morphological correlatio n between [Fe l l] and both the 
[O III] and radio 2-cm emission. iFalcke et al. Nl996ll present, 
in addition, an Ha map which shows a similar morphology. 
The integrated [Fell] fiux within a 5" radius is 4.43 x 10~^* 
erg s~^ cm~^. 

Colour map (J-K): The J ~ K colour map of the inner 
10" X 10" region of the galaxy is shown in Fig.|21 The main 
feature is a colour gradient from J ~ K — 0.95 to the SW to 
J — K = 1.34 to the NE. J — K ^ 1 is the color of a late-type 
stellar population characteristic of galaxy bulges. Assuming 
that the redder color is due to obscuration by dust, we esti- 
mate an average reddening Ay ~ 1.2 mag for the NE side. 
This result also suggests that the NE is the near side of the 
galaxy, where dust in the galaxy disk obscures the light from 
the bulge. In addition, there is a red feature extending from 
the nucleus to the N, possibly due to a dust lane. We esti- 
mate Av « 1.8 mag for this feature. The orientation of the 
photometric major axis obtained from the J-band image is 
129 °± 2. 
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Figure 2. Comparison between infrared, optical line emission and radio continuum maps. The solid straight line is the PA of the [Fell] 
maximum extent, (a) [Fell] (thin contours and gray scale) overlaid with the infrared J band image (thick gray contours). The flux is 
given in log(erg cm~^ pix~^) according to the gray scale bar at the top. (b) [O III] contours and gray scale in arbitrary units, (c) 
Radio 2-cm contours and gray scale in arbitrary units. The black circle indicates the region affected by seeing and the box represent the 
GNIRS IFU field of view. The images were rotated to coincide with the orientation of the IFU. 





5 



(arcsec) 

Figure 3. J — K colour image (gray scale and thin contours). 
The relation between J — K colour and shade of gray is given by 
the bar at the top. The J-band contours are shown overlaid in 
thick gray lines. The black circle indicates the region affected by 
seeing and the box represent the GNIRS IFU field of view. The 
images were rotated to coincide with the orientation of the IFU. 



3.2 2D intensity maps 

The final IFU data cube contains 520 spectra, each cover- 
ing 0'.'l5x O'.'lS, corresponding to 14 pc^ at the galaxy. The 
total IFU field of 3'.'9x3'.'0 thus corresponds to a region of 
dimensions 360pcx275pc at the galaxy. We fitted Gaussian 
profiles to the [FeiijA 1.257, H2A 2.121, Pa/3 and Br7 emis- 
sion lines in order to obtain the integrated fiuxes, radial 
velocities (from the peak wavelengths of the lines) and ve- 
locity dispersions (from the widths of the lines). A A 2.12/im 



continuum image was obtained from the interpolation of the 
continuum under the H2A 2.121 emission line. 

In Fig.0]we present the 2D maps of the near-IR con- 
tinuum fiux and emission line intensities in arbitrary units, 
with mean uncertainties of 5, 12, 15 and 17% for H2A 2.121, 
Br7, Pa/3 and [Fe 11] A 1.257, respectively. In the same figure 
we also present the H2/Br7, Br7/Pa/3 and [Fen]/Pa/3 line 
ratios. The nucleus in these maps is defined as the peak 
of the continuum light distribution and is identified by the 
white cross in each panel of Fig.^] In order to investigate 
the relation between the radio and near-IR emission, we have 
overlaid contours of the radio 2 cm continuum emission from 
iFalcke. Wilson fc Simpson 1 (Il99^ on the Br7 and [Fe 11] in- 
tensity maps as well as on the line ratio maps. In order to 
align the radio image with our images, we have adopted as 
the nucleus in the radio map the peak just to the left of the 
hot spot, as in Fig. H which we then registere d to the peak 
of our continuum image. This was also done bv lFalcke et al. I 
(jl99^ for registering their optical and radio images. They 
claim that the uncertainty between the relative position of 
the nucleus in the optical and radio bands is «0'.'2. In our 
case, the uncertainty may be smaller, due to the smaller ef- 
fect of reddening in the near-IR relative to the optical, but 
0'.'2 corresponds to only one pixel of our data, thus we may 
consider such small difference as a coincidence when com- 
paring our images with the radio images. 

The H2 intensity could be measured over most of the 
IFU field, while there was not enough signal to measure 
the Br7 and Pa/3 in a few lines close to the top and bottom 
edges ( l'.'5 from the nucleus). The [Fell] emission- line could 
only be measured over the inner l'.'9 along the minor axis 
of the IFU field, and thus the corresponding map covers a 
smaller region than those in the other emission lines. All 
emission lines are most extended along PA« 129°, which is 
the orientation of the radio jet (and galaxy line of nodes), 
and present a bipolar structure extended to both sides of 
the nucleus, in good agreement with our narrow-band [Fe 11] 
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H2A2.121/i.m Flux Continuum Image 
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Figure 4. Left panels show, from top to bottom, the intensity maps for H2A2.121^m, Br7, Pa/3 and [Fe ll] A 1.257/im emission lines 
with mean uncertainties of 5, 12, 15 and 17%, respectively. The right panels show, from top to bottom, the reconstructed image for 
the continuum emission at 2.12^m, H2A2.121/Br7 line ratio map, Br7/Pa/3 ratio map and [FeliJA 1.257/Pa/3 ratio map. The thin black 
contours overlaid on the H2 and Pa/3 intensity maps and on the continuum image are isointensity contours. T he thick white lines overlaid 
on the Br7, [Fell], H2/Br7 and [Fell] /Pa/3 maps are the contours of the 2cm radio continuum emission from lFalcke. Wilson fc Simpsoni 
^^9^ while on the Br7/Pa/3 map the contours are from the J- ii" color map of Fig.|3] The alignment uncertainty between the radio 
emission and our 2D maps is ~0''2. The spatial orientation and scale are the same for all figures and the cross marks the nucleus (the 
peak of the continuum emission) . 
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image, and with both the [O ill] and radio continuum images 
presented in Fig.|5| The overlaid contours of the radio con- 
tinuum image on the Br7 and [Fe ll] emission-line maps show 
a detailed correspondence between the radio and emission- 
line structures: the strongest line-emission is observed to the 
NW, and approximately coincides with the strongest emis- 
sion in radio; the emission-line distribution to the SE bends 
to the NE, as also observed in the radio map. The peak of 
line emission is nevertheless a bit displaced from the radio 
emission peak: the line emission peaks at O'.'GNW while the 
radio emission peaks at CK'S NW from the nucleus. 

A comparison between the emission-line maps of the H I 
with that of [Fell] shows that the latter has the strongest 
emission displaced to the NW relative to the Hi emission, 
tracing the structure of the radio hot spot which bends to 
the W. 

In the case of the H2 emission map, although the above 
bipolar structure is also observed, there is additional emis- 
sion at the lowest intensity levels extended beyond the bipo- 
lar structure. 



3.3 Line ratio maps 

In order to investigate the main excitation mechanisms of 
H2 and [Fell], we have constructed the line-ratio maps 
[Fell] /Pa/3 and Il2/Br7 shown in Fig.^] The uncertainties 
on these ratios range from 0.2 at the nucleus to 0.7 at the 
top and bottom borders of the IFU field for the first ratio 
and from 0.1 to 0.9 for the second. The uncertainties are 
thus quite large at the top and bottom regions of the line 
ratio maps. We have overlaid the radio contours also on the 
line ratio maps. Along the radio structure, the ratio H2/Br7 
is approximately constant, with a value 0.8 ±0.1. Outside 
the radio emission region, the H2/Br7 ratio increases, reach- 
ing values of 5.3 ± 0.9 at l'.'2 from the nucleus in the direc- 
tion perpedicular to the radio axis. The ratio [Fell]/Pa/3 
is 0.9 ±0.3 at the nucleus, and increase up to ^1.5 at the 
borders of the radio structure, except to the South of the 
nucleus, where it reaches the lowest value of 0.5 ± 0.2. From 
the emission-line maps, it can be observed that this region 
presents faint levels of both [Fe 11] and radio emission. 

We have constructed also the line ratio map Br7/Pa/3, 
which can be used as a reddening indicator. As we do not 
have a relative calibration between the J and K bands, this 
map can only be used to investigate the relative reddening 
distribution. The lowest ratios are observed to the SW, while 
the highest ratios are observed to the NE, the transition 
occurring approximately at the major axis of the galaxy. 
This behaviour is the same as the one observed in the J — K 
map of Fig.|3 suggesting that the line emission is subject 
to the same extinction as the continuum, with the NE more 
reddened than the SW side, consistent with the NE being 
the near side of the gala^xy. In order to evidence this, we have 
overlaid the J — K contours of Fig.|21on the Br7/Pa/3 ratio 
map, showing that the reddest J — K values correspond to 
the highest Br7/Pa/3 ratios. 



3.4 Gaseous kinematics 

Gaseous radial velocity fields were obtained measuring the 
peak wavelength of the emission lines using Gaussian curves 



to fit the emission lines. The results are shown in the left 
column of Fig.|K| where the red colors represent redshifts and 
blue colors represent blueshifts. The mean uncertainty in the 
velocities for all lines is less than lOkms"^. An underlying 
"rotation pattern" is present in the four maps although it is 
quite clear that in all cases there are other important kine- 
matic components, evidenced by large deviations from sim- 
ple rotation. We have overlaid the radio contours on the Br7 
and [Fell] radial velocity maps, which evidence the strong 
influence of the radio jet on the gas kinematics. In partic- 
ular, the location of the radio hot spot corresponds to the 
regions of highest blueshift observed in the gas towards the 
NW, while there is also some correspondence between the 
radio emission and the redshifts to the SE, mainly observed 
in the Hi and [Fell] velocity maps. 

The velocity field most closely resembling circular rota- 
tion (the classical "spider diagram" ) is the one derived from 
the H2 line. We have tentatively modeled the H2 velocity 
field as disk rot ation in a central Plummer bulge potential 
JPlummeJll91lll : 



GM 



CI 



(1) 



where M is the mass, Co the scale length, r is the radius 
in the plane of the galaxy and G is the Newton's gravita- 
tional constant. Defining the coordinates of the kinematical 
center of the system as {Xo,Yo), the observed radial veloc- 
ity at the position {R,"^), where R is the projected radial 
distance from the nucleus in the plane of the sky and ^ is 
the corresponding position angle, is given by 



Vr = Vs + 



/GMR sinicos^/^ ; 



{i?2[sin2(* - *o) + C0S2 i COs2(* - I-q)] ± COs2 



(2) 



where 14 is the systemic velocity, i is the inclination of the 
disk (j = for face on disk) and ^0 is the position angle 
of the line of nodes iStorchi-Bergmann. Wilson fc BaldwinI 
1996fl . The equation above contains seven free parameters, 
including the kinematical center, that can be determined by 
fitting the model to the observations. This was done using a 
non linear least-squares algorithm, in which initial guesses 
are given for the free parameters. The isovelocity curves rep- 
resenting the best model fit to the data are shown as black 
lines in the upper left panel of Fig.|^ while the white con- 
tours trace the isovelocity curves of the observed velocity 
field. It can be seen that circular rotation is not a good rep- 
resentation of the velocity field, particularly along the radio 
axis, evidentiating the presence of other kinematic compo- 
nents in the gas besides rotation. The kinematic center ob- 
tained from the fit coincides with the location of the peak 
of the continuum light distribution within the uncertainties. 

The parameters derived from this fit are: the systemic 
velocity - 1752 ± 2.8 km s~^ is 54kms~^ larger than the 
one listed in NASA/IPAC Extragalactic Database (NED), 
while the inclination of the disk - 54.8° ± 3.9° is 10° higher 
than the photometric one we obtain from the J-band image 
(cos~^ (^), where a and 6 are the semi-major and semi-minor 
axes of the ellipse fitted to the outermost isophote of Fig.^J. 
We obtain the value 1.02±0.13 x 10® Mq for the bulge mass 
and 108.5 ± 12.9 pc for the bulge scale length. The orienta- 
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tion of the line of nodes - 119.5° ± 1.4° is 10° smaller than 
that of the photometric major axis 

We have also measured the full width at half maximum 
(FWHM) of the observed emission lines, from which we ob- 
tained the velocity dispersion as ct = ^-y^H^. The corre- 
sponding a maps for each line are shown in the right panels 
of Fig.|S| The mean uncertainties in a values are 2, 4.5, 14 
and 17 kms~^ for H2, Br7, Pa(3 and [Fen], respectively. We 
have again overlaid the radio contours on the Br7 and [Fe 11] 
crmaps. It can be observed that the regions with the largest 
a are displaced relative to the regions with strongest radio 
emission: for example, the largest a values are observed be- 
tween the radio nucleus and the hot spot in the Br7 cr map. 
Similarly, in the [Fell] amap the highest values avoid also 
the locations of the radio emission knots. 

A comparison between the a maps of the different lines 
show that the H2 velocity dispersion values are lower than 
those of [Fell] and Hi lines, indicating that the gas emitting 
the two latter lines is more perturbed than the molecular 
gas. 

3.5 Emission-line "Tomography" 

The relatively high spectral resolution of the data has al- 
lowed us to slice our data cube into a sequence of velocity 
bins along each emission line profile, providing a better sam- 
pling of the gas kinematics, not restricted to the peak wave- 
length of the emission-line, but including also the wings. 
The slices were obtained after the subtraction of the contin- 
uum determined as averages of the fiuxes from both sides 
of each emission- line. Each slice corresponds to a velocity 
bin of «50kms~^, selected to include two pixels of the data 
cube and to correspond approximately to the nominal spec- 
tral resolution of the data. The obtained slices are shown in 
Fig.|S| where each panel presents fiux contour levels in arbi- 
trary units for each velocity slice and each emission line, and 
where the nucleus, defined as the peak of continuum emis- 
sion, is marked with a cross. The zero velocity is adopted as 
the one corresponding to the peak wavelength of H2 within 
an aperture of 3x3 pixels centered at the nucleus. We have 
overlaid the contours of the radio image on the central panels 
the Br7 and [Fell] emission line slices in Fig. |n] 

The velocity distributions are somewhat different for 
the different emission lines. For Hi, as the slices trace gas 
from negative (blueshift) to positive (redshift) residual ve- 
locities, the peak in the flux distribution moves from NW 
to SE, approximately following the path traced by the radio 
emission. The blueshifts observed in Hi, reach ~ 300 kms^^, 
approximately at the location of the radio hot spot, between 
0'.'6 and O'.'SNW of the nucleus. Some blueshifts are also ob- 
served to the SE following the radio structure which bends 
upwards (see Fig. |SJ. Redshifts are mostly observed at the 
nucleus, where they reach ~200kms~^ and to the SE where 
they reach smaller values of ~100kms~^ at 1-2" SE. The 
redshifted emission to the SE also traces the radio structure. 
Thus, to the SE, we observe both blueshifts and redshifts 
following the radio structure. 

For H2, the behaviour is similar to that observed for 
Hi, only that there is additional emission beyond the ra- 
dio structure. The highest blueshifts, of up to ~400kms~^ 
are observed in the [Fell] emitting gas, which shows a more 
bipolar structure in the velocity slices than the other emis- 



sion lines. Throughout this structure, we observe mostly 
blueshifts, even at the nucleus and to SE. 

In summary, we could say that the [Fell] velocity field 
shows predominantly a bipolar structure, which is also ob- 
served in the H2 and Hi for the highest velocity gas. For 
velocities between -100 and 100 kms~^, the structure in Hi 
is best described as linear, with both blueshifts and redshifts 
observed to either side of the nucleus. In the case of H2, in 
this low velocity range, the emission at low intensity levels 
is spread over most of the IFU field. 



4 DISCUSSION 
4.1 Gas kinematics 

Although all the velocity fields show evidence of a rota- 
tion pattern, the blueshifted side shows a more abrupt gra- 
dient going from the systemic velocity at the center to 
-lOOkms"^ at only 0'.'6 (55 pc) from it towards the NW, 
while to the SE equivalent redshift is barely reached at the 
border of the IFU field (210 pc from the nucleus). In the top 
panel of Fig.Qwe present one-dimensional cuts of the radial 
velocity field along the radio jet for the H2, Br7 and [Fell] 
emission lines, in which we can clearly see this behaviour. 

By comparing the line intensity maps (Fig.|lJ with the 
velocity fields (Fig.|5J , we verify that the peak blueshifts ap- 
proximately coincide with the peak intensities, suggesting 
that the flux enhancement is produced by the compression 
provi ded by the NW radio jet (Falckc, Wilson & Simpsoni 
Il99al . which should be at least partially oriented toward us 
in order to account for the observed blueshift. A more de- 
tailed analysis of the correspondence between the line fluxes 
and velocities (but along one axis) can be done using one- 
dimensional cuts shown in Fig.|7| by comparing the top and 
middle panels. In this comparison it can be observed that 
the [Fell] line flux presents a peak at the location of the 
strongest blueshift, between 0'.'6 and O'.'S NW of the nu- 
cleus, another smaller peak at the nucleus and another one 
at l'.'2 SE of the nucleus, at a location coincident with a 
"residual blueshift" - a decrease in the velocity observed in 
the redshifted side in the top panel of Fig.Q The fiuxes of 
the H2 and Br 7 emission lines also show peaks at «0'.'4 NW, 
thus somewhat closer to the nucleus than the locations of 
strongest blueshifts and between 1 and 2" SE of the nucleus, 
approximately at the locations of the residual blueshift de- 
scribed above, also observed in the H2 and Br 7 emission 
lines. Similar displacements be tween the radio and optical 
emission have been observed bv lFalcke. Wilson fc Simpson I 
I19981') in their comparison between HST [Olll]A5007 and 
[Nlls]-|-Ha emission- line images and the radio images of 
ESO 428-G14: they found that the emission-line fluxes are 
enhanced, in the form of knots and filaments in structures 
shaped as cocoons surrounding the radio jets. This close 
association between emission line distributions and radio 
structures is nevertheless not always observed leading some 
authors to conclude that there is no connection between the 
kinematics of the NLR and radio jets (Kaiser ct al. 200^; 
iDas et al E005. 2006). 

In the case of ESO 428-G14 the good correspondence be- 
tween the radio and emission line structures, also observed 
in the velocity slices of Fig. |S| indicates that the radio jet 
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Figure 5. Left: Prom top to bottom, velocity fields of tlie emitting gas for H2, Br7, Pa/3 and [Fell]. The mean uncertainty in velocity 
is less than 10 kms~^ for all lines. Right: Velocity dispersions maps for each emission line with mean uncertainties of 2, 4.5, 14 and 17 
kms~^ for H2, Br7, Pa/3 and [Pell], respectively. The thin white contours on the H2 and Pa/3 maps are isovelocity contours, while the 
heavy white contours show overlaid radio emission on the Br7 and [Pell] panels. The black lines on the H2 velocity map represent the 
modeling of the H2 velocity field by a Plummer potential. The spatial orientation and scale are the same as for Pig.BI 
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Figure 6. Velocity slices across the emission lines in ~50kms~^ velocity bins centered on the velocities indicated in each panel. The 
black cross marks the position of the nucleus (peak of continuum emission). Top: H2 (left), Br7 (right); Bottom: Pa/3 (left), [Fen] (right). 
We present flux contours for each slice except for the middle panels of Br7 and [Fell] where we present the radio continuum contours. 
The spatial orientation and scale are the same as for Fig. HI 



has an influence on all the emitting gas, although the some- 
what distinct kinematics observed in the different emission 
lines suggest varying contribution of the radio jet to the 
line emission. The [Fell] line kinematics is dominated by 
two outflowing structures towards the end of the radio jets; 
the H I kinematics has contribution from all emitting regions 
along the radio jet; and the H2 hue presents the less bipolar- 
ity, more closely resembling a rotation pattern. We interpret 
these differences as due to a larger disk rotation component 
contribution to kinematics of the H2 hue, an increasing con- 
tribution of the radio jet perturbations to the kinematics of 
the Hi line emission, and a even larger contribution from 
the radio jet to the kinematics of the [Fe 11] line. 

The influence of the radio jet on the emitting gas can 



also be observed in the velocity dispersion (a) maps (Fig.|^ 
and in the one-dimensional cuts along the radio axis shown 
in the lower panel of Fig[7| Increase in the a values are 
observed from 0'.'4SE to O'.'GNW along the radio axis and 
between O'.'S and I'.'S to SE of the nucleus. At least to the 
NW, the largest a values are thus observed between the 
nucleus and the radio hot spot, and do not coincide with 
the emission line flux and velocity peaks, but is shifted to- 
wards the nucleus. To the SE, the highest a values seem to 
approximately coincide with the peaks in the line intensity 
and velocity (the residual blueshifts relative to the redshifts 
of the gas to SE). We interpret this result as due to the in- 
teraction of the radio jet with ambient gas which produces a 
broadening of the emission lines due the momentum trans- 
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Figure 7. One-dimensional cuts along the radio axis. Top: Radial 
velocity for H2 Br7 and [Fell]. Middle: Emission line fluxes in 
arbitrary units. Bottom: Velocity dispersion values. 



fered to the gas by the jet. The different behaviour of the a 
values to the NW and to the SE can be understood as due 
to a shock between the radio jet and a dense ISM to the NW 
which blocks the radio jet and stirs the gas not only at the 
shock front but also behind the shock, which is observed as 
a larger velocity dispersion. The presence of a shock front is 
evidenced by the the abrupt termination of the radio jet at 
the hot spot and the enhanced line emission at this location. 
The interaction of the radio jet with the ISM is milder to 
the SE, where the radio jet stirs the gas but is not blocked 
penetrating farther into the ISM. 

It can also be noted that the overall velocity dispersion 
values are lowest for H2 (30 to 85 km s~^), and highest for the 
[Fell] (50 to 180kms~^), indicating that the [Fell] traces the 
most disturbed gas, while II2 traces the less disturbed gas. 
This result is also consistent with those obtained from the 
velocity fields discussed above, supporting a different origin 
for the gas emitting the different lines. Distinct kinematics 
for Hi, H2 and [Fell] has been obse rved also in previous 
near-IR studies o f other AGN s fe.g. IWilman et. al]l2000l 
IStorchi-B crgmann et al.lll999^ . where the broadest [Fen] 
emission has been attributed to gas cloud interactions with 
the radio jet. 



4.2 H2 emission origin 

The II2 lines can be excited in two different ways: fluorescent 
excitation through absorption of soft-UV photons (912-1108 
A) in the Lyman and Werner bands ( Black & van Dish oe^ 
Il987r) and collisional excitation, namely inelastic collisions 
between molecules in a warm (T > 1000 K) gas. The 
heating necessary to allow collisional excitation may be 
provided by: shocks jHollenbach fc McKed [l989^ . X-rays 



jMalonev. HoUenbach fc TielensI Il996l) or dense photodis- 
sociation regions by UV photons. The method commonly 
used to differentiate collisional excitation from fluorescence 
is based on the flux ratios of H2 lines in the K band. 
Based on measurements of such ratios for a lar g e sam ple of 
galaxies, iRodri'guez-Ardila. Riffel fc Pastorizal (|200^ con- 
cluded that fluorescence i s not important for AGNs support - 
mg previous studies by VciUoux. Goodrich fc Hill lll997t) . 
iMoorwood fc Oliyal Sw90 ) and Fischer ei _al,, (.1987,') . 
^"" louiflei^^nil^ggf) ■ using HST, imaged a sample of 
10 Seyfert galaxies in H2 and detected extended emission 
in six of them. For three galaxies they found H2 emission 
in the extended narrow line region, coincident with [O ill] 
and Hq-|- [N ii] line emission. For these galaxies they dis- 
carded UV excitation as the dominant excitation process 
on the basis of the Ha to H2 ratio, and argued that slow 
shocks are the most likely process to explain these line ra- 
tios. They com piled spectroscopic observa t ions of mole cular 
hydrogen fromp< !oornnee f fc Israeil lll996l) , iRuid lll997f) and 
IVeifleux. Goodrich fc Hid lll99'<D and looked for correlations 
with radio 6-cm and hard X-ray flux. They found no corre- 
lation with X-rays, suggesting this is not the dominant H2 
excitation mechanism, and found a weak correlation with 
radio 6-cm, suggesting that no single mechanism is likely 
to be responsible for the molecular hydrogen excitation in 
Seyfert galaxies. 

In a paper aimed at explain i ng th e strong H2 emission 
of NGC 6240, iDraine fc Woodj lll990D argued that shock 
speeds higher than the H2 a values of ES0 428-G14 (50^ 
a ^80kms~^) dissociates the H2 molecule. We thus con- 
clude that faster shocks, although present, as indicated by 
the larger [Fe 11] a values, destroy the H2 molecule, and this 
is why we do not see broader H2 emission. iDraine fc Wood^ 
(11990|) propose that most of the H2 line emission originates 
in molecular gas with is heated by transient X-ray irradia- 
tio n. Hard X-rays from the AGN have been also proposed 
bv lWilman et. alJ ll2000tl and iBeflamv fc Tadhuiiteil ^00^) 
as the dominant excitation mechanism of H2 emission in 
Cygnus A. 

Our H2 flux distribution and kinematics suggest some 
association with the radio morphology, in particular at 
the radio hot spot, but there is also additional emission 
spread throughout the IFU field. The emission line ra- 
tio H2A2.121/Br7 is «0.8 along the radio jet increas- 
ing outwards, in the direction perpendicular to the ra- 
dio jet, to ~5. In Starbursts, where the main heating 
agent is UV radiation, H2A2.121/Br7 < 0.6, while for 
AGNs this ratio is larger (0.6 < H2A2.121/Br7 < 2) be- 
cause of the ad ditional H2 ex cited b y shocks or by X-rays 
from the AGN dRodn'guez-Ardfla. Riffel fc Pastorizal 12003: 
IStorchi-Bergmann et al .ll999^ . The value « 0.8 observed for 
ESQ 428-G14 along the radio jet is thus typical of AGNs. 

The higher H2A2.121/Br7 line ratio values observed 
outside the jet could be due to some X-rays from the AGN 
still reaching these regions while the UV-photons needed to 
excite the H I are in much less number producing only very 
faint Br7 emission. 

From the discussion above we conclude that the colli- 
sional excitation necessary to produce the H2 emission may 
be provided by the interaction of the radio jet with the cir- 
cumnuclear interstellar medium (hereafter ISM) or by heat- 
ing produced by the X-rays emitted by the AGN. The en- 
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hanced H2 a values in the hot spot region reveal that the 
jet is depositing kinetic energy on the gas. This energy can 
be estimated as follows. The instrumental broadening is 21 
kms~^. Subtracting this in quadrature, we obtain a mini- 
mum a for H2 of 34 kms~^ (away from the radio jet), and a 
maximum value at the hot spot region of 77 kms~^ . (Uncer- 
tainties in these values are ~ 10 kms~^.) Velocities of the or- 
der of the latter provide enough kinetic energy to excite the 
H2 molecule. Assuming that the H2 emitting regions with 
smallest a (blue region in the top right panel of Fig. |^ are 
solely excited by X-rays and that the X-ray contribution is 
the same as in the jet region, we can estimate the increase in 
kinetic energy provided by the radio jet as (77/34)^ = 5. Un- 
der the above assumptions, one could conclude that ~83% 
(5/6) of the H2 excitation is due to shocks at the hot spot 
regions (red regions in the top panel of Fig.jKJ. In the yellow 
regions SE of the nucleus, and the regions surrounding the 
radio jet region, the kinetic energy ratio is 2.7, thus the jet 
excitation still dominates (~73%) there. But these are only 
upper limits because the X-ray flux in the jet region may be 
larger. In a "photoionization scenario" by a "clumpy" torus 
atmosphere (eg. ^yans ct al. 1993; Konigl & Kartie 1994), 
the ionization state of the gas would be larger along the 
coUimation axis of the torus and would drop with increas- 
ing distance from this axis. This could explain why the H2 
emission is distributed over a wider area, in gas which is 
irradiated by an attenuated continuum. 

4.3 [Fell] emission origin 

A number of papers have presented reviews and de- 
tailed studies of the physical conditions of the [Fen] 
emitting gas in s tarburst and AGN galaxies . suc h as 
Mouri ct al. ('l990VMouri. Kawara fc Tanieuchil l|l993^ and 
RodriEucz-Ardila, Riffcl fc Pastoriza (2005). In these stud- 
ies it is argued that the [Fell]/Br7 ratio is controlled by 
the ratio between the volumes of partially to fully ionized 
gas regions, as [Fell] emission is excited in partially ion- 
ized gas. Such zones in AGN are created by power law pho- 
toionization (including X-rays emitted by the AGN) or shock 
heating. These two processes are discriminated by the elec- 
tron temperature of the [Fe 11] region: Te « 8000 K for pho- 
toionization and Te ~ 6000 K for shocks. A contribution 
from shocks pro duced by radio jets is supported by results 
such as those of lForbes fc WardI l)l993f) which have found a 
correlation between the [Fe ll] and 6-cm radio em ission in 
radio galaxies. Nevertheless ISimpson et alj lll996l) have ar- 
gued that the dominant excitation mechanism of [Fen] is 
photoionization with shock excitation accounting for only ~ 
20% of the [Fell] emission in AGN. 

The nature of the excitation mechanism 
can be investigated using emission line ratios. 
[RodriEucz-Ardila, Riffcl & Pastoriza ll2005ll show that 
[Fell]/Pa/3 is smaller than 0.6 for Starburst galaxies and 
higher than 2 for supernova remnants, for which shocks 
are the main excitation mechanism. The Seyfert galaxies 
have values in between 0.6 and 2 suggesting that this 
ratio can be used as a measure of the relative contri- 
bution of photoionization and shocks. For ESO 428-G14 
l±0.27<[Fell]/Pa/3 <2±0.73 in the region co-spatial with 
the radio jet, indicating a range of relative contribution of 
photoionization and shock excitation for the [Fe 11] line, with 



regions where this ratio is « 2 been dominated by shocks. 
In the region corresponding to the radio hot spot this ratio 
ranges from 1.5 to 2 supporting a large contribution from 
shocks. 

In order to estimate an upper limit on the contribution 
of the radio jet in the excitation of the [Fell], we use the 
same reasoning as in the case of H2 emission above, assum- 
ing that the lowest a observed for this line (dark blue regions 
in the bottom right panel of Fig.|3 is produced by X-ray ex- 
citation. We obtain a ratio between maximum and minimum 
kinetic energies of (180/60)^ = 9, thus 90% of [Fell] emitted 
in the hot spot region (red region in bottom right panel of 
Fig-El is excited by the radio jet under the above assump- 
tions. In other regions of the jet (yellow-green regions) where 
a ^ 130kms~^, the percentage of gas emission due to radio 
jet excitation falls to 80%. But these are also upper limits 
because as pointed above for H2, the X-ray emission may 
be stronger along the radio jet, providing a larger contri- 
bution for the [Fen] emission excitation than in the regions 
with the lowest a values. On the other hand, the stronger 
association of the [Fen] flux and kinematics with the radio 
structure supports a larger contribution of the radio jet to 
the [Fen] excitation relative to that of H2. 



5 SUMMARY AND CONCLUSIONS 

We have analyzed near-IR J and K band 2D spectra of the 
inner Ri300pc of the Seyfert galaxy ES0 428-G14 obtained 
with the Gemini GNIRS IFU. The fine angular sampling 
of this instrument and the high spectral resolution of the 
observations provided a spatial sampling of 14 pc'^ at the 
galaxy and a velocity resolution of «20kms~^. These char- 
acteristics and 2D coverage has allowed us to obtain gas flux 
distributions and ratios as well as to map the gas kinemat- 
ics in the four strong emission lines [Fe n]A1.257/im, Pa/3, 
H2A2.121/im and Br7. Such 2D mapping at the spatial and 
spectral resolution of our observations and in the near-IR is 
unique in the literature so far. 

We have used the spectroscopic data in conjunction 
with a [Fe n]A1.257/xm narrow-band image and J — K color 
map obtained with OSIRIS at the CTIO Blanco Telescope. 
The [Fen] narrow-band image shows a bi-polar structure 
which correlates well with the structure observed in a pre- 
vious published radio map, as well as with that observed in 
a narrow-band [O ni] HST image. This bi-polar structure is 
oriented approximately along the photometric major axis of 
the galaxy at PA~129°. Under the assumption that redder 
color is due to dust screening, we conclude from the J — K 
color map that NE is the near side of the galaxy. 

The IFU observations show that not only the [Fe n] flux 
distribution, but also Pa/3 and Br7 correlate with the radio 
emission, all showing elongated morphologies along the radio 
axis. The H2 emission has an additional component more 
spread over the disk of the host galaxy. 

The gas kinematics shows a pattern which, at first sight, 
seems to be due to rotation in a disk, but which is clearly 
disturbed by non-circular motions. The analysis of the gas 
kinematics is complicated by the fact that the galaxy ma- 
jor axis is approximately aligned with the radio axis. The 
radio emission morphology suggests that a radio jet is the 
cause of outfiows observed in velocity "slices" obtained along 
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the gas emission line profiles: blueshifts of ~ 300-400 km s^^ 
are observed in emission associated with a radio hot spot at 
~0'.'8 NW of the nucleus, while mostly redshifts and some 
blueshifts observed to the opposite side of the nucleus are 
associated to fainter radio emission observed up to 2" SE. 
From this velocity pattern, as well as from the emission-line 
ratios we conclude that the radio jet leaves the nucleus at 
a small angle relative to the plane of the galaxy, with the 
NW side slightly directed toward us and the counter-jet di- 
rected away from us, being partially hidden by the disk of 
the galaxy. Some redshift observed to the NW and some 
blueshift observed to the SE shows that the line emission 
is observed to both sides of the plane of the galaxy, con- 
firming that the angle between the radio jet and the plane 
of the galaxy is small and emission associated with the jets 
occurs on both sides of the galaxy plane. The reddening dis- 
tribution obtained from the Br7/Pa/3 emission-line ratio co- 
incides with that obtained from the continuum, confirming 
that the line emission comes from circumnuclear gas located 
essentially in the galactic plane. 

The velocity dispersion maps show the highest values 
between the nucleus and the regions of strongest radio emis- 
sion (the radio hot spots). We interpret this result as due 
to kinetic energy deposited by ratio jet in the circumnu- 
clear ISM, producing a compression in the gas and emission 
enhancement just beyond this compressed region. 

While the [Fell] emission is dominated by the highest 
velocities at the two hot spots in opposite sides of the nu- 
cleus, the Hi emission has an important contribution also 
from lower velocities, observed all along the radio emission 
structure. The H2 emission has significant emission contribu- 
tion at low velocities spread throughout the observed region, 
apparently coming from the disk of the host galaxy. 

From the observed kinematics, we conclude that the 
radio jet has a fundamental role in shaping the emission 
line region as it interacts with the galaxy ISM surrounding 
the galaxy nucleus. We have used the 2D velocity dispersion 
maps to estimate the kinetic energy deposited in the circum- 
nuclear ISM by the radio jet relative to regions away from 
the jet. Assuming that the H2 and [Fell] excitation in the 
latter regions is dominated by X-rays, and that the X-rays 
excitation is the same in the jet region, we obtain contribu- 
tions by shocks of up to 80 - 90% for the [Fe 11] and up to 
70-80% for H2 in the jet region. These are however, upper 
limits due to the fact that the X-rays contribution along the 
jet axis may be larger than away from it. The stronger asso- 
ciation of the [Fe 11] emission and kinematics with the radio 
structure supports a larger contribution of the radio jet to 
the excitation of [Fell] than to that of H2. 
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